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Abstract 

The r~ — > r)ir~ TT + ir~ v T decay with the r\ — » 77 mode is studied using 384 ftr 1 of data collected 
by the BABAR detector. The branching fraction is measured to be (1.60 ± 0.05 ± 0.11) x 10 -4 . It 
is found that r _ — > /i(1285)7r~!/ T — > r;7r _ 7r + 7r _ i/ T is the dominant decay mode with a branching 
fraction of (1.11 ± 0.06 ± 0.05) x 10~ 4 . The first error on the branching fractions is statistical and 
the second systematic. In addition, a 90% confidence level upper limit on the branching fraction 
of the t~ — > ?/(958)7r~^ T decay is measured to be 7.2 x 10~ 6 . This last decay proceeds through a 
second-class current and is expected to be forbidden in the limit of isospin symmetry. 

PACS numbers: 13.35.Dx, 14.60. Fg 



I. INTRODUCTION 



The high-statistics sample of r-pair events collected 
by the BABAR experiment allows detailed studies of r- 
lcpton decays with small branching fractions. Many 
of these decays are poorly understood and more pre- 
cise measurements of the branching fractions as well 
as studies of the decay mechanisms are required. This 
work examines the r~ — > 7]n~ '-k + it~ v t decay [1] where 
7] — > 77. We show that this mode is dominated by 
t~ — > /i(1285)7r~^ T . This decay mode has been previ- 
ously studied by the CLEO collaboration [2] . A measure- 
ment of the t~ — + /i(1285)7t~j/ t decay was also made by 
the BABAR collaboration, but with a different final state 
[3]. 

This work also presents a search for the t~ — > 
rj' (958)ir~ v T decay where ?y'(958) — > i]Tr~ir + . Since this r 
decay proceeds via a second-class current, it is expected 
to be forbidden in the case of isospin symmetry. A 90% 
confidence level upper limit has been previously set by 
the CLEO collaboration at 7.4 x 10" 5 [2]. 

This analysis is based on data recorded by the BABAR 
detector at the PEP-II asymmetric-energy e + e~ storage 
rings operated at the Stanford Linear Accelerator Center. 
The data sample corresponds to an integrated luminosity 
of 384 fb _1 recorded at center-of-mass energies of 10.58 
GeV and 10.54 GeV between 1999 and 2006. With a 
cross section for e + e~ — > t + t~ production of (0.919 ± 
0.003) nb [4], this data sample contains approximately 
706 million r decays. 

The BABAR detector is described in detail in Ref. [5]. 
Charged particle momenta are measured with a five-layer 
double-sided silicon vertex tracker and a 40-layer drift 
chamber inside a 1.5-T superconducting solenoidal mag- 
net. A detector of internally reflected Cherenkov light 
(DIRC) provides tt/K separation. A calorimeter consist- 
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ing of Csl (Tl) crystals is used to measure the energy of 
electromagnetic showers, and an instrumented magnetic 
flux return (IFR) is used to identify muons. 

Monte Carlo simulation is used to evaluate the back- 
ground contamination and selection efficiencies. The 
r-pair production is simulated with the KK2F Monte 
Carlo event generator [6]. The r decays, continuum qq 
events, and final-state radiative effects are modeled with 
Tauola [7], JETSET [8], and Photos [9], respectively. The 
generic r Monte Carlo sample contains one decay mode 
with an 77 meson in the final state, r~ — > r]7T~ it v T . In 
addition, a dedicated Monte Carlo sample is generated 
using KK2F and Tauola for the r~ — > rjK®ir~v T decay. 

Dedicated samples of t + t~ events are created using 
EvtGen [11] where one of the r leptons can decay to 
any mode included in Tauola [7] and the other r de- 
cays to an r]ir~ ir + n~ v T final state. One of the samples is 
generated using the t~ — > /i(1285)7t _ j/ t decay, and the 
other is generated using t~ — > 7777" tt + tt~ v t phase space. 
The /i(1285) meson decay modes that are relevant to 
this analysis are the /i(1285) — > rjTr~ir + (40%) and the 
/i(1285)-> a (980)7r (60%) decay modes where the rel- 
ative contributions to the r~ — > /i(1285)7r _ ^ T decay are 
indicated in parentheses [10]. The Monte Carlo distri- 
butions identified as signal in the figures use a combina- 
tion of r~ — > i]iT~ 'n + ~tt~ v T phase space, and the resonant 
t~ — > /i(1285)7r _ ^ T decay samples where the relative 
fraction is based on the branching fractions measured in 
this work. 

The detector response is simulated with GEANT4 [12]. 
All Monte Carlo simulation events are passed through 
a full simulation of the BABAR detector and are recon- 
structed in the same way as the data [5]. 



II. SELECTION 

Events of interest are isolated with a loose pre- 
selection. Since r pairs are produced back-to-back in 
the e + e~ ccntcr-of-mass frame, the event is divided into 
hemispheres according to the thrust axis [13], calculated 
using all reconstructed charged and neutral particles. 
The analysis procedure selects events with one track in 
one hemisphere (tag hemisphere) and three tracks in the 
other hemisphere (signal hemisphere). The total event 
charge is required to be zero. 
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Charged particles are required to have transverse mo- 
menta greater than O.lGeV/c in the laboratory frame. 
The distance of the point of closest approach of the track 
to the beam axis must be less than 1.5 cm. In addition, 
the z coordinate (along the beam axis) of the point of 
closest approach of the track must be within 10 cm of 
the z coordinate of the production point. Neutral clus- 
ters are required to have an energy of at least 30 MeV 
and must not be associated with a charged track. 

After pre-selection a more discriminating analysis se- 
lection is applied. This selection strategy has three com- 
ponents. The first selection criterion is based on the 
event-shape properties. The magnitude of the thrust is 
required to be between 0.92 and 0.99, in order to re- 
duce the non-r backgrounds. The background from non- 
r sources can arise from Bhabha, di-muon, two-photon 
and qq events. 

The second set of selection criteria requires the parti- 
cles in the tag hemisphere to originate from a leptonic r 
decay (r _ — > e~ vv or t~ — > /i _ z/77). The track in the tag 
hemisphere must be identified as an electron or muon and 
must have a momentum in the center-of-mass frame be- 
low 4 GeV/c. The first criterion removes qq events while 
the second removes lepton-pair events. Electrons are 
identified with the use of the ratio of energy measured 
by the calorimeter to track momentum (E jp\ the ioniza- 
tion loss in the tracking system (dE/dx), and the shape 
of the shower in the calorimeter. Muons are identified 
by hits in the IFR and energy deposits in the calorime- 
ter expected for a minimum-ionizing particle. Residual 
background from qq events is reduced by requiring that 
there be at most one electromagnetic calorimeter cluster 
in the tag hemisphere with energy above 50 MeV and 
that the total neutral energy in the tag hemisphere be 
less than 1 GeV. 

The final set of selection criteria is applied to the sig- 
nal hemisphere. The aim is to reduce the residual back- 
grounds from r decays while maintaining high selection 
efficiency for t~ — > rjn~ tt + tt~ v t decays. An event is re- 
jected if any of the tracks in the signal hemisphere is iden- 
tified as an electron or if any pair of oppositely charged 
tracks is consistent with originating from a photon con- 
version. 

The event selection requires that there be one unique 
r\— > 77 candidate in the signal hemisphere. The 77—^ 77 
candidates consist of two neutral clusters in the electro- 
magnetic calorimeter with an invariant mass (M~~) be- 
tween 0.47 and 0.62 GeV/c 2 . To reduce combinatoric 
background from other r decays with ir° mesons, the 
higher and lower-energy clusters must have E > 0.7 GeV 
and E > 0.3 GeV, respectively. 

Residual background from other r decays and qq events 
is reduced by requiring that there be no 7r° mesons in the 
signal hemisphere, where a tt° candidate consists of two 
neutral clusters in the electromagnetic calorimeter with 
an invariant 77 mass between 115 and 150 MeV/c 2 . In 
addition the invariant mass of the 7/7r~7r + 7r~ system is 
required to be less than 1.8 GeV/ c 2 . No particle iden- 



tification algorithm is applied to the charged tracks to 
distinguish pions from kaons, and all invariant masses 
are calculated assuming that the charged particles are 
pions. Note that the t~ — > fi(1285)K~i> T decay is kine- 
matically disfavored. Monte Carlo studies indicate that 
the background from r decays is dominated by 77 com- 
binatorics in the 7r~7r + 7r~(> l7r°) mode. The absolute 
amount of r background is determined in the fits. The 
backgrounds from r~ — > riK2.ir~i> T and qq are evaluated 
separately and discussed in section III. 




Mass (GeV/c 2 ) 



FIG. 1: The invariant masses of the (a) 77, and (b) r)-K~TY + -K~ 
final states are shown. The dark shaded histograms show the 
simulated signal events, the lightly shaded histograms show 
the simulated r background and the unshaded histograms 
show the simulated qq background. All selection criteria are 
applied. In (b) the cut requirement on the invariant mass of 
the r]-K~TT + TT~ system is not imposed and the invariant mass 
of the 77 system is between 0.50 and 0.58 GeV/c 2 . 



III. RESULTS 

The invariant mass of the 7/ candidates is shown in 
Fig. 1(a) after all selection criteria are applied. The plot 
shows a difference in the reconstructed mass and width 
of the 77 meson between the data and Monte Carlo sim- 
ulation. The width of the 77 peak obtained in a fit using 
only the signal Monte Carlo sample was slightly narrower 
but consistent with the width from the data. The vari- 
ations in the fit results between data and signal Monte 
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FIG. 2: The invariant masses of the (a) r/n~TT + and (b) rjn^ 
final states are shown. The dark shaded histograms show the 
simulated signal events, the lightly shaded histograms show 
the simulated r background and the unshaded histograms 
show the simulated qq background. Note that (a) and (b) 
have two entries per event. All selection criteria are applied. 
In (b), it is required that the invariant mass of the associated 
t]ty~tv + system is between 1.23 and 1.32 GeV/c 2 . 



Carlo are accounted for with a systematic error discussed 
in the following section. 

Figure 1(b) shows the invariant mass of the ?77r~~7r + 7r~ 
system after all selection criteria (except on this variable) 
have been applied. The disagreement between the data 
and Monte Carlo simulation in Fig. 1(b) shows that the 
underlying physics is more complex than the model used 
to simulate the decay and may involve additional reso- 
nances. For example the model used in Tauola for the 
r — > 37r27r ^ T mode, does not give an accurate repre- 
sentation of the experimental data [14]. However, the 
?/7r~7r + 7r~ invariant mass distribution is not used in the 
determination of the branching fractions presented in this 
paper and the modeling uncertainties are small and in- 
cluded in the systematic errors. 

Figure 2(a) shows the invariant mass of the rjir~ir + 
system. Figure 2(b) shows the invariant mass of the 
7y7r~ and t]tt + systems with the requirement that the in- 
variant mass of the 777r~7r + system is between 1.23 and 
1.32 GeV/c 2 (consistent with being an /i(1285) meson). 
Only the ir mesons forming the /i(1285) candidates are 
shown in this plot. The peak at 980 MeV/c 2 is due to 



the a (980) in the /i(1285) -> a (980)7r decay. 



A. Inclusive r — > t]tt tt + it v t branching fraction 

The invariant mass distribution of the 77 system is 
fitted with a Novosibirsk function [15] (Gaussian distri- 
bution with a tail parameter) for the r\ meson plus a 
polynomial function for the background. The fit range is 
0.47 to 0.62 GeV/c 2 and the fit is a binned X 2 nt. The 
observed width of the 77 is dominated by the experimen- 
tal resolution (14 MeV/c 2 ). The resolution and the tail 
parameters in the fit to the data are fixed to the values 
obtained from a fit to the signal Monte Carlo simulation 
(see the following paragraphs for a discussion of the sys- 
tematic error associated with this constraint). The peak 
position and normalization parameters of the Novosibirsk 
function are allowed to vary in the fit to minimize the 
dependence of the result on the difference in the rj mass 
observed between the data and Monte Carlo simulation. 
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FIG. 3: Plot (a) shows the invariant mass of the 77 system. 
Plot (b) shows the invariant mass of the data rjn~n + system. 
The points represent the data, the solid line is the fit function, 
and the dashed line shows the background. 

A total of 2174 ± 73 events are obtained from the fit 
shown in Fig. 3(a). The x 2 / n -d-f for the fit is 51/46. 

The t~ — > 777T - 7r + 7r~ v T branching fraction is measured 
with 
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where N b s is the number of events obtained from the 
fit, Nbkgd is the number of background events with an 77 
meson (371 ± 83), N T + T - is the number of r leptons in 
the sample calculated from the luminosity and e+e~ — > 
t+t~ cross section, e is the efficiency for selecting the 
signal events (4.18 ± 0.06%), and B(rj^ 77) is 0.3943 ± 
0.0026 [10]. The t~ — > -qn^ n + ir~ v T branching fraction 
is measured to be (1.60 ± 0.05 ± 0.11) x 10~ 4 where the 
first error is statistical and the second is systematic. 

The systematic errors are dominated by the uncer- 
tainty in the results of the fit to the 77 meson (5.0%), 
which is partly due to the difference in the 77 mass res- 
olutions in the data and Monte Carlo. The sensitivity of 
the results to the fit to the 77 peak is investigated by un- 
constraining the width and tail parameters in the fit to 
the data. Also, polynomials of different orders are tested 
as background functions. The 5.0% uncertainty associ- 
ated with the fit is also due partly to the variation of 
the branching fraction observed for different background 
functions. The remaining systematic errors include terms 
for the uncertainties of the 77 background levels (3.8%), 77 
selection efficiency (3.0%), track reconstruction (2.4%), 
lepton identification (1.6%), selection efficiency statis- 
tical error (1.4%), luminosity (1.0%), and the 77 — > 77 
branching fraction (0.7%). 

The background events {Nbkgd) includes a contribution 
from qq events which is estimated from qq Monte Carlo 
samples. The uncertainty in the number of background 
events extracted from the qq background is evaluated by 
comparing data and Monte Carlo simulation distribu- 
tions in regions where there is an enhanced amount of 
qq events (events with an ■q-n~'K + 'K~ invariant mass that 
is larger than the r mass). The qq Monte Carlo predicts 
125 candidates with an 77 meson with an uncertainty of 
18 events. 

The background also gets a contribution from t~ — > 
i]K®ir~is T events. The number of K® candidates is de- 
termined by counting the number of events that pass the 
full selection from the dedicated r~ — > rjK^-K^v T Monte 
Carlo. The Monte Carlo predicts 246 background 
events. The uncertainty on the number of selected back- 
ground events is dominated by the uncertainty of the 
t~ — > r]K a s r K^v T branching fraction B T -^ nK a 7T -^ — 

(1.10 ± 0.35 ± 0.11) x 10~ 4 [16]. The total background 
is estimated to be Nbkgd = 371 ± 83 where the uncer- 
tainties from the qq and K° backgrounds are added in 
quadrature and included as a systematic error. 

The stability of the branching fraction measurements 
was tested by varying the selection criteria (within a 
range of values determined by the level of agreement be- 
tween data and Monte Carlo), which did not change the 
results significantly. Furthermore, the results of the fit 
to the 77 meson mass peak are found to be stable to vari- 
ations in the bin width or mass range used in the fit. 
Branching fractions measured on the electron and muon 
samples separately are found to be consistent. 



B. Branching fractions for 
t~ — > /1 (1285)7T _ v T — > »77r _ 7r + 7r - v T and 
T" -> / 1 (1285) 7 r-iy T 

The /i(1285)7r candidates are selected with the crite- 
ria specified in the previous section with the requirement 
that the 77 candidates have 0.50 < M 77 < 0.58 GeV/c 2 . 
The invariant mass of the rjTT~ir + system is fitted with a 
Breit-Wigner function convoluted with a Gaussian distri- 
bution and the background is modeled with a Novosibirsk 
function summed with a polynomial (see Fig. 3 (b)). The 
X 2 /n.d.f. is 180/152 for the 777r"7r+ fit. A P-wave Breit- 
Wigner function [17] is used to fit the data while the fit 
to the simulated distribution uses a simple Breit-Wigner 
function as implemented in the generator. In both cases 
the Breit-Wigner function is modulated by phase space. 
The normalization and mean of the Breit-Wigner func- 
tion are allowed to vary and the width is fixed to the 
PDG value. The resolution parameter of the Gaussian 
function is fixed to 7 MeV/c 2 , which is the mass resolu- 
tion of the r]ir~n + system obtained from simulation. 

The background function is determined by fitting the 
777r~7r + invariant mass distribution obtained from a sam- 
ple of simulated t~ — » r/ir^ ir + n~ v T events where the de- 
cay does not proceed through an /i(1285) meson. 

There is no evidence for the production of /i(1285) 
mesons in the data from background sources. This is 
determined by relaxing selection criteria to increase the 
background from multihadron events. 

The r~ — » /1 (1285)71-^1/,- — > r/7r~ n + n~ v T branching 
fraction is determined with 

Nobs 11 
Br-^n-^-u^ fl ) = 2N T+T -~eB{7 1 ~, li y (2) 

where N D b s is the number of /i(1285) mesons obtained 
in the fit (1255 ±70), N T + T - is the number of r pairs ob- 
tained from the luminosity and e + e~ — ► t + t~ cross sec- 
tion, e is the efficiency for selecting ar^ /i(1285)7r~^ T 
event (4.08 ± 0.07)%, and B(r\ —> 77) is the 77 -> 77 
branching fraction (0.3943 ± 0.0026) [10]. The r" -> 
/i(1285)7r~z^ r branching fraction is determined by divid- 
ing Eq. 2 by the fi(1285) -> 777^+ rate (0.35 ± 0.11) 
[10]. 

The branching fractions for the t — » /i(1285)7r v t — > 
777r~ 'ir + ir~ v T and t~ — > / 1 (1285)7r _ ^ T modes arc (1.11 ± 
0.06 ±0.05) x 10~ 4 and (3.19 ± 0.18 ± 0.16 ± 0.99) x 10~ 4 
respectively, where the first error is statistical and the 
second is systematic. The third error quoted on the 
t~ — ► /i(1285)7T _ ^ T measurement is due to the large er- 
ror on the /i(1285) — > i]-k^tt + branching fraction. Most 
systematic errors for these branching fractions are com- 
mon to the ones listed for the inclusive measurement. 
While the 77 fit uncertainty affects the inclusive result 
only, an extra systematic error of 1% comes through 
the /i(1285) decay modeling due to the uncertainty of 
the branching fractions of the /i(1285) — > a (980)7r and 
/i(1285)^ 777r~7r + decay modes [10]. This is determined 
by varying the relative contribution of the two modes 
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within the quoted uncertainties. 

The fraction of the t~ — > /i(1285)7r~i/ T — > 
?77r~7r + 7r~^ T mode to the inclusive t~ — > rjTT~ tt + tt~ v T 
mode is found to be 0.69 ± 0.01 ± 0.05 where the first 
error is statistical and the second is systematic (taking 
into account the correlations between the various com- 
ponents). 

C. Limit on the r~ — > r]' (958)tt _ v-t branching 
fraction 

A limit on the r~ — ► 7/(958)7r _ kv branching fraction 
can be set by searching for decays of the rf (958) to the 
r)Tr~Tr + final state. This r decay mode proceeds through 
a forbidden second-class current and is not expected to 
produce an observable signal [18]. 
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FIG. 4: Fit to the r/(958) region of the r/-K^TT + invariant 
mass spectrum. The fit uses a Gaussian distribution for the 
peak summed with a linear function for the background. The 
events are tagged with muons or electrons, and all selection 
criteria have been applied. 

A fit to the ?77r~7r + mass distribution is performed with 
a Gaussian function for the rf (958) and a polynomial 
function for the background (see Fig. 4). The mean of 
the Gaussian is fixed to the mass of the 7/(958) meson. 
The width of the Gaussian distribution is fixed to the 
value obtained in a fit to a data sample containing a 
significant number of 77' (958) mesons. This data sample 
is created by removing all selection criteria except the 
loose pre-sclcction described in section II. 

We observe 19 ± 13 candidates. To set a limit, we 
treat all of the events in the 7/ (958) peak as signal. We 
assume that the efficiency for selecting t~ — > r]' (958)w~ v T 
events is the same as the r~ — ► /i(1285)7r _ z^ T selection 
efficiency. The systematic uncertainty is dominated by a 
7% error due to the uncertainty in the mass resolution 
of the 7/ (958). The remaining systematic errors are the 
same as those described in the previous section. The 
results give a 90% confidence level upper limit on the 
T ~ - * rf (958)n~ v T branching fraction of 7.2 x 10~ 6 . 



IV. SUMMARY 

The r~ — > rjn~ n + n~ v T decay using the 77^ 77 mode 
is studied with the BABAR detector. It is found that 
r~ — > fi(1285)ir~ v T is the dominant decay mode for the 
r]Tr~7t + iT~ final state. 

The branching fraction of r~ — ► r/7r~7r + 7r~^ T is mea- 
sured to be (1.60 ± 0.05 ± 0.11) x 10~ 4 where the first 
error is statistical and the second systematic. This mea- 
surement is more precise than the CLEO result (2.3 ± 
0.5) x 10~ 4 [19]. 

The branching fraction of the r~ — ► /i(1285)7T _ z>y — > 
?77r~ 7r + 7r~ v T decay mode is measured to be (1.11±0.06± 
0.05) x 10~ 4 and is consistent with previous results [10]. 

The branching fraction of t~ — > /i(1285)7r~i/ r is mea- 
sured to be (3.19 ± 0.18 ± 0.16 ± 0.99) x 10~ 4 . The 
first error is statistical, the second is systematic, and the 
third error is associated with the 30% uncertainty on the 
fi — > 7/7r~7r + branching fraction [10]. This measurement 
is in agreement with the CLEO result of 5. 8^3 x 10~ 4 [2] 
and the BABAR result of (3.9 ± 0.7 ± 0.5) x 10~ 4 [3]. The 
branching fraction of r~ — > fi(12&h)ir~v T is predicted by 
effective chiral theory to be 2.9 x 10~ 4 [20]. 

A 90% confidence level upper limit on the branching 
fraction of the r~ — > rf' (958)n~ v T decay is measured to 
be 7.2 x 10~ 6 . This is an order of magnitude lower than 
the previous 90% confidence level upper limit of 7.4 x 
10~ 5 set by the CLEO collaboration [2]. No significant 
evidence for this second-class current decay mode of the 
t is observed. 
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